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Abstract 

In the present work we considered Galilean conformal algebras 
(GC A) , which arises as a contraction relativistic conformal algebras {xi 
exi, t ^ t, e — > 0). We can use the Galilean conformal symmetry to 
constrain two-point and three-point functions. Correlation functions 
in space-time without boundary condition were found in [191 . In real 
situations there are boundary conditions in space-time, so we have cal- 
culated correlation functions for Galilean confrormal invariant fields in 
semi-infinite space with boundary condition in r = 0. We have calcu- 
lated two-point and three-point functions with boundary condition in 
fixed time. 
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1 Introduction 



Recently, there has been some interest in extending the AdS/CFT corre- 
spondence to non-relativistic field theories [U [2]. The Kaluza-Klein type 
framework for non-relativistic symmetries, used in Refs. [Hill [3], is basically 
identical to the one introduced in [1] (see also [5]). The study of a different 
non-relativistic limit was initiated in [6] , where the non-relativistic conformal 
symmetry was obtained by a parametric contraction of the relativistic con- 
formal group. Galilean conformal algebra (GCA) arises as a contraction rel- 
ativistic conformal algebras [71 E] ; where in d = 4 the Galilean conformal 
group is a fifteen parameter group which contains the ten parameter Galilean 
subgroup. Beside Galilean conformal algebra, there is another Galilean al- 
gebra, the twelve parameter schrodinger algebra [H [2] . The dilatation gen- 
erator in the schrodinger group scales space and time differently, Xj — )■ Axj, 
t — )■ X'^t, but in contrast the corresponding generator in GCA scales space 
and time in the same way, Xi — > Xxi, t ^ Xt. Infinite dimensional Galilean 
conformal group has been reported in [iTj , the generators of this group are : 
= -(n + l)t''xidi - r+^dt, = f'+^di and J^- = -f^i^Xidj - xjdi) for 
an arbitrary integer n, where i and j are specified by the spatial directions. 
There is a finite dimensional subgroup of the infinite dimensional Galilean 
conformal group which is generated by (J?-, L^^, L'^, M^^ , M^). These gen- 
erators are obtained by contraction ( t — )• i, Xj — )■ exj, e — >• 0, ~ e ) of 
the relativistic conformal generators. Several conformal extensions of the 
Galilean Lie algebra have been obtained recently 0. 

The representation for finite GCA was found in [lOJ and pjj. According 
to the theorem [12] "Every representation D{G) of a finite group on an 
inner product space is equivalent to a unitary representation," the repre- 
sentations of finite Galilean conformal group are unitary and we can use 
this group for physical applications. For example, non-relativistic limit of 
conformal hydrodynamic describes the small fluctuations from thermal equi- 
librium. Whenever the conformal invariant dissipative hydrodynamics are 
considered in nonrelativistic limit, the relations d^T^^'^ = and = ( 
The stress-energy tensor of the CFT obeys these relations ) are converted 
to incompressible Navier-Stokes (NS) equation [13]. From reference [H] 
it can be inferred that, the Navier-stokes equation for incompressible flow 
(V.v = 0, V is velocity of fluid ) is covarianted under the infinite GCA. The 
NS equation for ideal hydrodynamics (without viscosity) changes to Euler 
equation. The finite Galilean conformal group is symmetry of the Euler 
equation [7] (see also [151 [16]). The gravity dual of finite GCA was consid- 
ered in [71[H1[IZ] and the metric with finite 2d GCA isometry was obtained in 
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|18j . We can use the Galilean conformal symmetry to constrain correlation 
function [19]. Correlation functions of Galilean conformal invariant fields 
without any boundary condition were found in [19]. 

The presence of free surfaces or walls in macroscopic systems which are at 
the critical point, lead to the large variety of physical effects. Since, using 
boundary condition effects is shown to be very helpful in various branch in 
physics, the systems with boundary conditions have been considered by both 
theorists [20] and experimentalists [21\ . The situation with walls or free sur- 
faces opens a new area in condensed matter physics [22]. In reference [23], 
the research on semi-infinite systems which exhibits a non-equilibrium bulk 
phase transitions was initiated and the effects of boundary condition on di- 
rect percolation were considered. 

Correlation functions near the boundary are different from other places, and 
in real situation there are boundary conditions in space coordinates. By 
using some methods in non-relativistic conformal field theory, we have ob- 
tained an analytic expression for correlation functions with boundary condi- 
tion. The situation when a system is in a predefined initial state and relaxes 
toward its critical equilibrium considered as a situation with a boundary 
condition at fixed time [23], for these reasons, in this paper we calculated 
two-point and three-point functions in space-time with boundary condition. 
The discussion will be exclusively in two dimensions, but the extension to 
arbitrary dimension is immediate. A system with boundary condition at a 
surface r = was kept invariant under the transformations were generated 
by [L_i, Lq, Li], but space translations, Galilean Boost transformation and 
space special Galilean conformal transformation no longer leave the system 
invariant. A system with boundary condition at a fixed time was kept invari- 
ant by the subalgebra [M_i, Mq, Mi, Lq]. Correlation functions with these 
boundary conditions are different with the corresponding results found from 
Galilean conformal invariance without boundary condition [19]. If the do- 
main of space-time coordinates are infinite in extent (without boundary con- 
dition) and the scalar fields are invariant under the Galilean conformal (GC) 
group, the two-point and three-point functions are completely determined 
[l9] . For example, these results apply to the calculation of time-delayed 
correlation functions of non-relativistic (small viscosity) systems at equilib- 
rium and at static critical point. If the space geometry is semi- infinite, a 
form of correlation function will be derived, see Eqs. (fT2]) and (p7|) . These 
results may be relevant to critical dynamics close to a surface, see |24j for 
an example, when a system is in a predefined initial state, it can be seen 
that critical relaxation towards equilibrium displays scaling at intermediate 
time p3j. (This situation is discribed as a system with boundary condition 
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in time coordinate.) We calculate the form of correlation functions with this 
condition, see Eq. (pOj) . 

The paper organized as follows. Section 2 is a brief review of GCA. In 
section 3, we calculate the form of correlation functions with a boundary 
condition at surface r = . Then in section 4, we calculate the form of 
two-point function in space-time with a boundary condition at a fixed time. 
Then in section 5, we extend these calculation to the three-point correlation 
function. Finally, in section 5, we close by some concluding remarks. 



2 Representations of Galilean conformal group 

Galilean conformal algebras (GCA) was obtained via a direct contraction of 
conformal generators. Physically, this comes from taking t ^ t , Xi ^ exi 
where e — >■ 0. The generators of conformal group are 

P^ = -id^, D = -ix^df,, J^J.u = i{x^du - x^d^), (1) 
= -iilx^x'd, - x^d^) 

(where fj,,^ = 0, 1, d). 
From the above scaling we obtain the Galilean conformal vector field gener- 
ators 

L_i = dr M_i = dt Lo = {tdt + rdr) Jij = -{xidj - Xjdi) (2) 
Li = {2trdr + t^dt) Mi = -t^dr Mq = -tdr 

Af_i, Lq) Jij: Mq, Li, Ml, are spatial translations, time translation, 
dilatation, rotations, boosts, time component and space components of spe- 
cial conformal transformation respectively. In following, the discussion will 
be exclusively in two dimensions, but the extension to arbitrary dimension is 
immediate. The above generators obey the following commutation relations 
in two dimensions, where define the Galilean conformal algebras. 

(m - n)Lm+n (3) 
{m - n)Mm+n 




[Lrm Ln] — 
[Lm,Mn] = 
[Mn,M.ra\ = 
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The above symmetries use to constrain two-point function (see section (3, 
4)). The representations of this group was obtained in |10j . According to 
the theorem [12| "every representation D(G) of a finite group on an inner 
product space is equivalent to a unitary representation" , the representations 
of 2d finite Gahlean conformal group are unitary and it can be used for 
physical applications. The representations of Galilean conformal group are 
built by Hausdorff formula [19] . 

[L^i,^]=dt<P (4) 

[Lo,(j)] = {tdt + rdr + A)(l) 

[Li,(j)] = {t^dt + 2trdr + 2tA - r^)(j) 

and 

[M_i,(A] = a,0 (5) 

[Mo,4>] = {-tdr + 04> 
[Mi,4>] = {-t^dr + 2tO^ 

where A is scaling dimension and ^ is rapidity . 

3 Two-point function in semi-infinite space 

Galilean conformal symmetry is related to massless non-relativistic sys- 
tems. Correlations of Galilean-quasiprimary fields were studied without any 
boundary condition [19], but in real systems with boundary conditions, we 
are interested in correlation function near the boundary, so we consider the 
effect of surface at r = 0. It is kept invariant under the transformations were 
generated by [L^i, Lq, Li], but space translations, Galilean Boost transfor- 
mation and space special Galilean conformal transformation no longer leave 
the system invariant. Nevertheless, it is known that Galilean conformal 
invariance can be used in analogous situations to constrain the two-point 
function |19] . For Galilean-quasiprimary fields, we require covariance only 
under the subalgebra [L_i,Loi^i]- 

Consider the two-point function of Galilean-quasiprimary fields 

G = G{ra,n,ta,h) =< 4>a{ra,ta)4>b{rh,th) > (6) 

and we require space points to be in the right half-plane, i.e. ra,rfe > 0. 
Time translation invariance gives G = G{ra,ri),T), with t = ta — tj,. From 
scale invariance we obtain 



5 



<0\[Lo,(t>a4>b]\0>=0 (7) 
i=b 

i=a 

= {Tdr + r„a,„ + ndr, + A)G = 

where A = Aq + A5. On the other hand, from the invariance under the time 
special GaUlean conformal transformation we have 



<0I |0>=0 (8) 

i=b 

^ ^{iidu + 2tiridi + 2tiAi - riii)G 

i=a 

= {{tl-tl)dr + 2{taradr, + hnOr,) 
+2{taAa + tbAh - Vaia + nib))G 
= {r^dr + 2tb{Tdr + Vadr^ + r^^rj 

+2{raCa - nCa) + ^TVadr, + 2(t^A„, + tbAb))G 

= {T^dr + 2{raia " TfoCa) + 2rr«a,, + 2tAJG = 

where in the last equation the scale invariance of G was used. 
Now, we make the ansatz 



G{ra,H,T) = T-^^^G'{u,v), tx = ^, v='^ (9) 

r T 

which solves for scale invariance, while Eq. ([8]) gives 

{udu - vd, - 2{uia + vib))G'{u, v) = 0, (10) 

Note that, the change of variables in Eq.(l9]) are not singular, because r = 
ia — tb 7^ 0. The general solution of Eq. (|lU|) is found by the method of 
characteristics ([20] and [25]). 

G'iu, v) = x{uv) exp(2(nCa - vi^)) (11) 
where x is ^in arbitrary function. The final result is 

G(r„r,,T) = 5A,,A,T-2^"X(^)exp(-(r,ea-rb6)) (12) 
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It is clear that, two-point function near the boundary is different from other 
places. We note that analogously to the conformal result |19j . the scaling 
dimension have to agree, while in this case we do not have a constrain on 
the rapidity ^a, since the system is not space special Galilean conformal 
invariant. One can use the relation (|12p for nonrelativistic conformal hy- 
drodynamics with small viscosity near the boundary. Two-point function 
in the bulk (out of boundary) for nonrelativistic conformal hydrodynamics 
was found in [19j . 



4 Two-point function for a non-stationary state 

We now consider a situation with a boundary condition at a fixed time. 
Boundary conditions of this type are kept invariant by the subalgebra [M_i, Mq, 
Ml, Lq\. For example, this may correspond to the situation when a system is 
in a predefined initial state and relaxes toward its critical equilibrium state 
|23j . Consider two-point function of quasiprimary fields 

G = G{ra,n, ta, h) =< <i)a{ra, ta)<i)b{n, U) > (13) 

Invariance under space translation implies G = G{r,ta,tb) with 
r = ra — Vf). We next demand invariance under Galilean Boost transforma- 
tion. 

<0| [Mo,./.a.^d |0>=0 (14) 

i=b 

Y^i-tidr, + Ci)G 

i=a 

= {-tadva +ia- tbdr^ + S,b)G 
= {-{ta-h)dr+^a + ^b)G = 

Space special Galilean conformal invariance demands that 

<0\[Mi,(Pa(pb]\0>=0 (15) 

i=b 

^{-t^idr, + 2Uii)G = 

i=a 

which gives respectively 

= 6 (16) 
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So, the two-point function reads: 

G = G'{ta,h)eM^^) (17) 

ta — tb 

In analogy to what was done before, we demand invariance under scale 
invariance 

<0l [Lo,0a06] |0>=0 (18) 
i=b 

i=a 

= {tadu + radu + Aa + tfe^t, + ^dt^, + ^b)G 

= {tadu + hdt, + rdr + Aa + Afe)G = 

we find 

G'(t„,t,)=t-(^^+^^)cI>(^) (19) 

tb 

where ^ is an arbitrary function. The final result is 

G{ta, t,) = t-^'^^^'^Mr) exp(-^) (20) 

tb ta — tb 

Note that here we have no condition on the exponents because the system 
is not invariant under the time special Galilean conformal transformation. 
We are able to use the above relation for nonrelativistic conformal hydrody- 
namics with small viscosity in critical equilibrium [23]. 
Our study in sections 3, 4 was an extension of the work of Henkel in [20] . 
This paper has devoted to the study of the schrodinger symmetry algebra, 
the maximal symmetry of free schrodinger equations, and one of the things 
has done there was considering two-point functions invariant under a sub- 
algebra of the schrodinger algebra. Here we applied the construction to the 
Galilean Conformal Algebra. Our results in Eq. (12) and Eq. (20) are cor- 
respond to the Eq. (3.35) and Eq. (3.44) respectively in [20]. The difference 
between these results is the form of the exponential piece ( r/t in the case 
of the GCA as opposed to r^/i in the schrodinger algebra) which is due to 
the difference in the form of the dilatation generator. 
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5 Three-point function in space-time with bound- 
ary condition 



Three-point function for Galilean conformal invariant fields in the bulk with- 
out boundary condition was calculated in [19]. In this section, we calculated 
three-point function near the boundary for Galilean conformal invariant 
fields. Three-point function in real situation with a boundary condition 
at fixed time was calculated. Along lines similar to section (3), we wish to 
construct three-point function of three GC-invariant fields 0j(Aj,,^j) with 
i = a,b,c. As we see in section (3), the correlation functions near the 
boundary are kept invariant under the transformations were generated by 
[L^i, Lq, Li]. Consider the three-point function of GC-invariant fields 

G{ra, n, rc, ta,tb, tc) =< (pa{ra,ta)(l)birb, tb)4>c{rc: tc) > (21) 

G is invariant under the time translation which is generated by Lq, so G = 
G{ra,rb,rc,T,a) where t = ta — and a = tb — t^. From scale invariance 
one can obtain 

<0| [Lo,0a0fe0c] |0>=0 (22) 

i=c 

^{tidt^ + Tidi + Ai)G 

i=a 

= {rdr + ada^ + Vadr^ + ndr^, + Tc^r, + + Afe + Ac)G = 

From the invariance under the temporal non-relativistic special conformal 
transformation we have 

<0I [Li,0„0fe0j |0>=0 (23) 

i=c 

^ ^(.tldu + '^ti^idi + 2UAi - riii)G 

i=a 

= {{tl - tc)dr + {tl - tl)d^ + 2{taradr^ + hndr, + tcTcdc) 
+2{taAa + tbAb + tcAc) - - H^b " rc^c)G 
= {r'^dr + a'^da- + 2tc{Tdr + Crda + radr^ + rfoSrj + rcdc) 
-raia - nia " T^ic + 2Tradr, + 2arbdr^ + 2{taAa + hAb + tcAc))G 
= {r'^dr - Taia " Uia " r^c + 2Tradr^ 

Wd^ + 2arbdr, + 2TAa + 2aAb)G = 

where in last equation the scale invariance ()22p of three-point function G 
was used. The above equations have this solution 

G' = .5A„+A„A,T-2^'"a-2^^G' (24) 
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If we introduce G' = G[{ra,rc,T)G2{rb,a) or G' = G[{n,rc,a)G2{ra,T), 
the above equations can be simplified as 

(T^dr + 2TTadr, - T^ia " r,i,)G'^ = (25) 
{Tdr+radr^+r,drJG[=0 

{a^ + 2arbdr, - n(b)G'2 = 

{ad^ + ndr,)G'2 = o 

or 

{a^d^ + 2andr, - nib - rcic)G'^ = (26) 
{ada + ndrf^ + rcdrJG'i = 
(t2 + 2rradr^ - raia)G'2 = 
(ra, + r,a,JG'2 = 

The general solution of these equations are found by using the method of 
characteristic |12j . 

G = 5A.+A„A.(ta - ^-'''^{tb - Q-'^^ exp(-^ + -^) (27) 

'■a T^c i^b ''c 
II ' a' c \ I ' d^c \ , I ' b' c \ I ' c^c \\ 

+T,{exchanging b c or a -H- c) 

where xi ^^'i X2 are arbitrary functions. Three-point function near the 
boundary is different from other places. In the following, we calculate three- 
point function for a situation with boundary condition at fixed time. Similar 
to section (4) three-point function (j2ip is kept invariant by the subalgebra 
[M-i, Mq, Ml, Lq]. Invariance under the spatial translation which is gener- 
ated by M_i implies G = G{r, s, ta,tb, tc) with r = — Tc and s = Vb — Vc- 
Invariance under non-relativistic Boost transformation is demanded 

<0| [Mo,(Aa</'fe</'c] |0>=0 (28) 

i=c 

^Y.^-Udr,+i^)G 
i=a 

= {-tadra +ia- hOr^ + 6 " ^c^r, + Cc)G 

= {-{ta - tc)dr - {tb - Qds +ia + ib + QG = 

SO, three-point function reads 

G = G'(t„4,t,)exp(-^ + '^k±^) (29) 
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Space special Galilean conformal invariance demands that 

<0\[Mi,(l)a(l)bM\0>=0 (30) 

i=c 
i=a 

which gives 

ic = 26 ia = ih (31) 

Finally we demand invariance under non-relativistic dilatation 

<0 I [Lc-^a'^foc/'J I 0>=0 (32) 

i=h 
i=a 

= {tadta + radra + + hdtt + Hdr-f, + Afc + tcdt, + rcdr, + Ac)G 

= {tadt, + tbdt, + tcdt, + rdr + sds + + A;, + Ac)G = 
{tadt^ + hdt, + t^dt, + A„ + Afe + A,)G' = 

we find 

G' = t-(^''+^'>+A.)cI>i(^) +t-(^'^+^<'+^=)cD2(^) +t-(Aa+A,+A.)^^(ia^^33^ 

where <l'i,<I>2j and $3 are arbitrary functions. The general solution is 
G[ra,n,rc,ta,tb,tc) = <^g.,gb%,g.+g6exp[ _ — + — — J 

[^-{Aa+A6+Ae)^^(^^) _^ ^^(Aa+At, + Ac)^^j-^^ _^ ^-(A, + A6 + Ae)^^|-^^j 

+T, (ex changing b ^ c or a ^ c) 

The above result is different from three-point function of the GCA without 
boundary condition. 

6 Conclusion 

We can use finite Galilean conformal group to constrain two-point and three- 
point functions. Correlation functions of Galilean conformal invariant fields 
in the bulk out of boundary were found in [19j . Correlation functions near 
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the boundary are different from other places, and in real situation there are 
boundary conditions in space coordinates. When a system is in an initial 
state and relaxes toward its critical equilibrium considered as a situation 
with a boundary condition at fixed time. In this paper we considered two 
real situations: 

1. A system with boundary in space coordinate (r = 0) was considered in 
sections (3), (5) 

2. A system with boundary condition at fixed time was considered in sec- 
tions (4), (5). 

The main results of this paper are the explicit expressions for Galilean con- 
formal invariant correlation functions in a semi-infinite geometry as given 
in Eqs. ([27]) and We calculated two-point function with 

boundary conditions at fixed time and surface r = 0, the form of two-point 
functions (jl2p . (j20p obviously are different with the corresponding results 
found from Galilean conformal invariance without boundary condition |l9j. 
The form of three-point functions are obviously different with corresponding 
results found from GC-invariant without boundary condition jl9j . 

Since Galilean conformal symmetry is related to massless (small viscos- 
ity) non-relativistic systems, one of the applications of GCA is considering 
nonrelativistic conformal hydrodynamics with small viscosity. One can use 
the Eqs. (jl2p . (j27p and (j20p . (j34p for nonrelativistic conformal hydrodynam- 
ics with boundary condition in surface r = and fixed time respectively. 
We know that in 2-dimensional space, there is a so-called exotic central 
extension of the GCA [26] . Martelli and Tachikawa [27] looked into the two- 
point functions in the bulk, but what happens close to a surface no-one has 
studied yet, as far as we know. We keep this interesting study for our future 
work in this topic. 

7 Acknowledgments 

We thank Prof. Malte Henkel for reading the paper and helpful comments 
and suggestions. 

References 

[1] D. T. Son, Phys. Rev. D 78, 046003 (2008). 

[2] K. Balasubramanian and J. McGreevy, Phys. Rev. Lett. 101, 
061601,(2008). 



12 



[3] C. Leiva and M. S. Plyushchay, Ann. Phys. 307, 372, (2003). 

[4] C. Duval, G. Burdet, H.P. Kiinzle, and M. Perrin, Phys. Rev. D31 1841 
(1985). 

[5] C. Duval, G. W. Gibbons, and P. A. Horvathy, Phy. Rev. D43, 3907 
(1991); C. Duval, M. Hassaine and P. A. Horvathy, Ann. Phys. 324, 
1158, (2009). 

[6] J. Lukierski, P. C. Stichel and W. J. Zakrzewski, Phys. Lett. A 357, 1 
(2006); J. Gomis, J. Gomis and K. Kamimura, JHEP 0512, 024 (2005). 

[7] A. Bagchi and R. Gopakumar, JHEP 0907, 037, (2009). 

[8] C. Duval and P. A. Horvathy, J. Phys. A42, 465206 (2009). 

[9] M. Henkel, Nucl. Phys. B 641, 405, (2002); M. Henkel, and J. Unter- 
berger Nucl. Phys. B 660, 407-435 (2003); R. Cherniha, M. Henkel, J. 
Math. Anal. Appl, 369, 120, (2010); M. Henkel, S. Stoimenov, Nucl. 
Phys. B847, 612, (2011); S. Fedoruk, P. Kosinski, J. Lukierski, and P. 
Maslanka (2011), [ar Xiv: 1012. 0480] : S. Fedoruk, E. Ivanov, J. Lukierski 
(2011), [a rXiv: 1101. 1658] , 

[10] M. Henkel, R. Schott, S. Stoimenov, and J. Unterberger, (2006), 
| arXiv:math-ph/06010 28| . 

[11] J. Negro, M. A. del Olmo and A. Rodriguez-Marco J. Math. 
Phys. 38, 3786, and 3810 (1997); C. Duval, P. Horvathy, (2011), 
jarXiv:1104.1502]. 

[12] Wu-Ki Tung, Page 35, Theorem 3.3, World Scientific, (1985). 

[13] I. Fouxon and Y. Oz, Phys. Rev. Lett. 101, 261602 (2008). 

[14] A. Mukhopadhyay, JHEP 1001, 100 (2010). 

[15] S. Bhattacharyya, S. Minwalla and S. R. Wadia, JHEP 0908, 059, 
(2009). 

[16] L Fouxon and Y. Oz, JHEP 0903, 120, (2009); P. M. Zhang, P.A. 
Horvathy Eur. Phys. J. C 65:607-614, (2010). 

[17] M. Alishahiha, A. Davody and A. Vahedi, JHEP 0908, 022, (2009). 

[18] A. Bagchi, A. Kundu, (2010), [arXiv: 101 1.4999) [hep-th]]. 



13 



[19] A. Bagchi and I. Mandal, Phys. Lett. B 675, 393-397, (2009). A. 
Bagchi, (2010), [ arXiv:1012.33l6 j: M. R. Setare, V. Kamali, (2010) 
^arXiv:1010.0329 [hep-th]]. 

[20] M. Henkel, J. Statist. Phys. 75 (1994). 

[21] K. Binder, P. C. Hohenberg, Phys. Rev. B6, 3461 (1972); T. C. Luben- 
sky, M. H. Rubin, Phys. Rev. Lett 31, 1469 (1973); T. C. Lubensky, M. 
H. Rubin, Phys. Rev. B 11, 4533 (1975); R. Lipowsky, Phys. Rev. Lett 
49, 1575 (1982); H. W. Diehl, C. Domb, J. L. Lebowitz, Vol. 10, p. 75. 
London, Academic Press (1986). 

[22] K. Binder, C. Domb, J. L. Lebowitz, Vol. 8, p. 1-144. London, Academic 
Press (1983) 

[23] H. K. Janssen, B. Schaub and B. Schmittmann, Z. Phys. B73, 539 
(1989). 

[24] M. Droz, L. Frachebourg and M. C. Marques, J. Phys. A24, 2869 (1991) 

[25] E. Kamke, Vol. 2, 4th edition, Akademische Verlagsgesellschaft (Leipzig 
1959). 

[26] C. Duval and P. A. Horvathy, Phys. Lett. B 479, 284 (2000); J. Lukier- 
ski, P. C. Stichel, and W. J. Zakrzewski, Phys. Lett. B 650, 203 (2007) 

[27] D. Martelh; Y. Tachikawa, JHEP 1005, 091, (2010). 



14 



